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The rare semileptonic B c — > D*t + 1~ decay is studied in the scenario of the universal extra dimen- 
sion model with a single extra dimension in which inverse of the compactification radius R is the only 
new parameter. The sensitivity of differential branching ratio, total branching ratio, polarization 
and forward-backward asymmetries of final state leptons, both for muon and tau, to the compacti- 
fication parameter is presented. For some physical observables uncertainty on the form factors and 
resonance contributions have been considered in the calculations. Obtained results, compared with 
the available data, show that there appear new contributions due to the extra dimension. 
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I. INTRODUCTION 



Flavor-changing neutral current (FCNC) b — > s,d transitions which occur at loop level in the standard 
model (SM) provide us a powerful tool to test the SM and also a frame to study physics beyond the 
SM. After the observation of b — > S7 these transitions became more attractive and since then rare 
radiative, leptonic and semileptonic decays of B u ^, s mesons have been intensively studied [2]. Among 
these decays, semileptonic decay channels are significant because of having relatively larger branching 
ratio. The experimental data for exclusive B — > K^*'£ + £~ also increased the interest in these decays. 
These studies will be even more complete if similar studies for B c , discovered by CDF Collaboration [3|, 
are also included. 

The B c meson is the lowest bound state of two heavy quarks, bottom b and charm c, with explicit flavor 
that can be compared with the cc and bb- bound state which have implicit flavor. The implicit-flavor states 
decay strongly and electromagnetically whereas the B c meson decays weakly. B UJ i, s are described very 
well in the framework of the heavy quark limit, which gives some relations between the form factors of the 
physical process. In case of B c meson, the heavy flavor and spin symmetries must be reconsidered because 
of heavy b and c. On the experimental side of the decay, for example, at LHC, 10 10 i? c events per year 
is estimated [3]- Hi- This reasonable number is stimulating the work on the B c phenomenology and this 
possibility will provide information on rare B c decays as well as CP violation and polarization asymmetries. 

In rare B meson decays, effects of the new physics may appear in two different manners, either through 
the new contributions to the Wilson coefficients existing in the SM or through the new structures in the 
effective Hamiltonian which are absent in the SM. 

Considering different models beyond the SM, extra dimensions are specially attractive because of in- 
cluding gravity and other interactions, giving hints on the hierarchy problem and a connection with string 
theory. Those with universal extra dimensions (UED) have special interest because all the SM particles 
propagate in extra dimensions, the compactification of which allows Kaluza-Klein (KK) partners of the 
SM fields in the four-dimensional theory and also KK modes without corresponding the SM partners ■ 
Throughout the UED, a simpler scenario with a single universal extra dimension is the Appelquist-Cheng- 
Dobrescu (ACD) model [l(|. The only additional free parameter with respect to the SM is the inverse of 
the compactification radius, 1/R. In particle spectrum of the ACD model, there are infinite towers of KK 
modes and the ordinary SM particles are presented in the zero mode. 

This only parameter have been attempted to put a theoretical or experimental restriction on it. Tevatron 
experiments put the bound 1/R > 300GeU. Analysis of the anomalous magnetic moment and B — » X s "f 
also lead to the bound 1[R> 300GeU. In the study of B — > K*-f decay [lH, the results restrict R to 
be 1/R> 250GeU. Also, in [13[ this bound is 1/R> 330 GeV. In two recent works, the theoretical study 
of B — >• AT 777 matches with experimental data if 1/R > 250 GeV [l4| and using the experimental result 
[l5j and theoretical prediction on the branching ratio of A& — > Afj, + fx~, the lower bound was obtained to 
be approximately 1/R ~ 250 GeV [Hj]. In this work, we will consider 1/R from 200 GeV up to 1000 GeV, 
however, under above consideration 1/R = 250 — 350 GeV region will be taken more common bound region. 
In literature, effective Hamiltonian of several FCNC processes [lj], Gil, semileptonic and radiative decays 
have been investigated in the ACD model [l9i - [29| . 



Concentrating on B c — > _D*l + l~decay, it has been studied by using model independent effective Hamil- 
tonian [3(|, in Supersymmetric models [3l| and with fourth generation effects [34]. Also in [33|, the UED 
effects on branching ratio and helicity fractions of the final state D* meson were calculated using the form 
factors obtained through the Ward identities for this process. The weak annihilation contribution in addi- 
tion to the FCNC transitions was taken into account. We will, however, only consider the FCNC transitions 
and calculate the lcpton asymmetries adding the resonance contributions. 

The main aim of this paper is to find the effects of the ACD model on some physical observables related 
to the B c — > D*£ + £~ decay, while doing this we also give the behavior of these observables by a couple of 
figures in the SM. Measurement of final state lepton polarizations is an useful way in searching new physics 
beyond the SM. Another tool is the study of forward-backward asymmetry (Apg), especially the position of 
zero value of Ap b is very sensitive to the new physics. In addition to differential decay rate and branching 
ratio, we study forward-backward asymmetry and polarization of final state leptons, including resonance 
contributions and uncertainty on form factors in as many as possible cases. We analyze these observables 
in terms of the compactification factor and the form factors. The form factors for B c — > Dl£ + £~have been 
calculated using the light front, constitute quark models l34j. the relativistic constituent quark model [35| . 
relativistic quark model [36| and light-cone quark model j37j. In this work, we will use the form factors 
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calculated in three-point QCD sum rules (38[. 

The paper is organized as follows. In Sec. II, we give the effective Hamiltonian for the quark level process 
b —> s£ + £~a.nd mention briefly the Wilson coefficients in the ACD model; a detailed discussion is given in 
Appendix A. We drive matrix element using the form factors and calculate the decay rate in Sec. III. In 
Sec. IV, we present the forward-backward asymmetry and Sec. V is devoted to lepton polarizations. In 
the last section, we introduce our conclusions. 



II. EFFECTIVE HAMILTONIAN AND WILSON COEFFICIENTS 



The quark-level transition of B c — > D*£ + £ decay is governed by b —> s£ + £ and given by the following 
effective Hamiltonian in the SM 39] 



V2n 



H eff = -^v tb v; a 



C% ff (s^L b) £Y l + C w (s ltl L b) £Y l ld 



(1) 



where q is the momentum transfer, L,R = (1 ± 7s)/2 and CiS are the Wilson coefficients evaluated at the 
b quark mass scale. 

The coefficient Cg^ has perturbative and resonance contributions. So, Cg can be written as 

Cf(/i)=C 9 M(l + ^( S ')) +Y(»,s') + Cr(»,s') (2) 

where s' = q 2 /m%. 

The perturbative part, coming from one-loop matrix elements of the four-quark operators, is 
Y(jjl,8') = h{y,s')[3C 1 (ti) + C 2 (fi) + 3C 3 (m) + C 4 (j«) + 3C 5 (m) + CfeO*)] 

- s') (4C 3 (/x) + 4C 4 (m) + 3C 5 ( M ) + C 6 (^)) 

- ~h(0,s') [C 3 (m)+3C 4 (m)] 

+ ~ (3C 3 (m) + C 4 ( M ) + 3C 5 (^) + C 6 ( M )) , (3) 
with y = m c /mb. The explicit forms of the functions w(s') and h(y, s') are given in [40|-[4"lj|. 

The resonance contribution due to the conversion of the real cc into lepton pair can be done by using a 
Breit-Wigner shape as [42j], 



ul m t ^ sml - m\r. + imvXv 



x [3d(/i) + C 2 Qi) + 3C 3 (/x) + C 4 (/x) + 3C 5 (^) + C 9 (jj)] . (4) 

The normalization is fixed by the data in 43] and the phenomenological parameter n is taken 2.3 to produce 
the correct branching ratio BR(B -» J/ijiK* -» K*£ + £~) = BR(B -» J/iiK*)B(J/tp -» £+£'). 

In the ACD model, there are not any new operators, therefore, new physics contributions appear by 
modifying the Wilson coefficients available in the SM. In this model, the Wilson coefficients can be written 
in terms of some periodic functions, as a function of compactification factor 1/R. The function F(xt, ^/R) 
which generalize the Fo(x t ) SM functions according to 

oo 

F(x t ,l/R) = F {x t ) + Y / F n(xt,x n ) (5) 

n=l 

where xt = mf/m^Y, x n — m^jm^ with the mass of KK particles m n — n/R. n=0 corresponding the 
ordinary SM particles. The modified Wilson coefficients in the ACD model, taken place in many works in 
literature, are discussed in Appendix A. 
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Briefly, for Cg, in the ACD model and in the NDR scheme we have 



C 9 ( M , l/R) = P NDR + Y{Xt, 2 l /R) - *Z(x t , l/R) + P E E(x u l/R). 
Instead of Cj, a normalization scheme independent effective coefficient Cy can be written as 

cfViAR) =v 16/23 c 7 (vw,i/R) 

a 8 
+ 1^14/23 _ r/ 16/23 )Cg(m; 1/R) + C2(pwi 1/R) J- htf*. 



The Wilson coefficient Cio is independent of scale fi and given by 

Y{x u l/R) 



C w (l/R) = 



i=i 



(0) 



(7) 



(8) 
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FIG. 1. The variation of Wilson coefficients with respect to l/R at q = lAGeV for the normalization scale 
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fi = 4.8 GeV. (Cg does not include resonance contributions. 



The Wilson coefficients differ considerably from the SM values for small R. The variation of modified 
Wilson coefficients with respect to l/R at q 2 = 14 GeV 2 , in which the normalization scale is fixed to 

fi = fib ~ 4.8 GeV, is given in Fig. [TJ The suppression of for l/R = 250 — 350 GeV" amount to 

75% - 86% relative to the SM value. |C 10 | is enhanced by 23% - 13%. The impact of the ACD on 
is very small. For 1/R> 600 GeV the difference is less than 5%. 



a 



eff 
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III. MATRIX ELEMENTS AND DECAY RATE 



The hadronic matrix elements in the exclusive B c — > D*£ + £~decay can be obtained by sandwiching the 
quark level operators in the effective Hamiltonian between the initial and the final state mesons. The 
nonvanishing matrix elements are parameterized in terms of form factors as follows (441 . |45| 



(D*(p Dt ,e) \s lfl (l - l5 )b\ B c (p Bc )) = -e^e* v p^- 



rn B + m D , 



is*(m Bc + m D ,)A 1 {q 2 ) 



and 



+* (PB e + PD t )„ (e*q) M f ] + iq, (e* q) ^ [A 3 (q 2 ) - A, (q 2 )} , (9) 
m Bc + m D * q z 



(D*(p D: ,e) 1^^(1+75)61 B c (p Bc )) = 2e^ aP e w p a D ,q fi T 1 {q 2 ) 



- m D ») - (p Bc +p D ')^(e*q) T 2 {q ) + i(e*q) q^ - (p Bc + p D *) 



/ M 2 9 

m Bc ~ m D* B 



Uq 2 ), (10) 



where q = p Ba — p B * is the momentum transfer and e is the polarization vector of D* meson. 
The relation between the form factors Ai(q 2 ), A 2 (q 2 ) and A 3 (q 2 ) can be stated as 



A (q 2 ) — mBs + m ^ A '" 2 ^ rnBs ~ m ^ 



Ivtia 



'-Mq 2 



Ivtia 



A 2 (q 2 



and in order to avoid kinematical singularity in the matrix element at q 2 — 0, it is assumed that 
4,(0) = MO) and Ti(0) = T 2 (0) @. 

Using the effective Hamiltonian and matrix elements in Eqs. (|9j)- (|10[) . the transition amplitude for 
B c — > D*£ + £~is written as 

M(B C -+ D* s £+£-) = ^V tb V* 

x J lfi\- 2Ae^ a ps w p a D ,qP - iBe*+iC(e*q)( PBc + PDj ) M + iD(e*q)q^ 



2Ee^ afj e w Pc ,,y - iF E * + iG{e*q){p Bc +Pd : )» + ^H{e*q)q, L 



(11) 



with the auxiliary functions 



A = C, 



eff 



V{q 2 ) , 2m b 



+ ^C^^T.iq 2 ), 



m Bc + m B * q 



B = C; ff ( mBc+mD *JA 1 ( q 2 ) + ^C e /H< c -™l*)T2(q 2 ), 

q 2 " 



c = a 



eff 



A 2 {q 2 ) , 2m b 



m Bc + m D , q 



: + ^C^(T 2 (q 2 ) 



n(q 2 )), 



D = 2C e g f ^(A 3 (q 2 ) A (q 2 )) 2^C^T 3 (q 2 ), 



q 

V(q 2 ) 



E — Ciq- 

m Bc + m D * 

F = C 10 {m Bc +m D - : )A 1 {q 2 ), 
A 2 (q 2 ) 



G = d 



m Bc + m D * 



H = 2C 10 ^(A 3 (q 2 )-A (q 2 )). 

q 



(12) 



Integrating over the angular dependence of the double differential decay rate, following dilepton mass 
spectrum is obtained 

dT G 2 a 2 m Ba . *, 2 rr A / 1Q s 



where s = q 2 /m Br , A = 1 + r 2 + s 2 — 2r — 2s — 2rs, r = m 2 D ,jm 2 B ^ v = Jl — 4m 2 /sm B ^ and 



A ^ K 

A D » =-Am Bc s 
2 



(3 - « 2 ) |A| 2 + 2w 2 |£| 2 



l\ m 2 Bc (3-v 2 )\C\ 2 +m 2 BcS 2 (l-v 2 )\H\ 2 



(3 - v 2 )(r + s - 1) - 3s(l - v 2 ) i?e[FG*] + 2m 2 Bc s(l - r)(l - u 2 )i?e[GF*] 



-2s(l - « 2 )i?e[FiP] + -(3 - w 2 )(r + s - l)Re[BC*\ 



Am|j A-3s(s-2r-2)(l- v 2 ) |G| 2 + A + 24rsv 2 |F' 2 



(A + 12rs)|B| z 



(14) 
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FIG. 2. (color online) The dependence of differential branching ratio on s, including the uncertainities on form 
factors in non-resonance case. (In the legend 1/R = 200, 350, 500 GeV.) 
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FIG. 3. (color online) The dependence of differential branching ratio on s with the central values of form factors 
including resonance contributions. 



In the numerical analysis, we have used ra^ = 6.28 GeV, mjj* = 2.112 GeV , mi, = 4.8 GeV, m M = 
0.105 GeV, m T = 1.77 GeV, \V tb V* s \ = 0.041, G F = 1.17 x 10~ 5 GeV^ 2 , r Bc = 0.46 x 10~ 12 s, and the 
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FIG. 4. (color online) The dependence of differential branching ratio on 1/R, with and without resonance contri- 
butions, including uncertainty on form factors at s = 0.18 for jj,, and s = 0.4 for r. (The subscript R in the legend 
represents resonance contributions.) 



values that are not given here are taken from [431 ] . In our work, we have used the numerical values of the 
form factors calculated in three point QCD sum rules [38| . in which q 2 dependencies of the form factors 
are given as 



F(q 2 



F(0) 



i + a( q y m % c ) + b( q ymij 2 ' 

and the values of parameters F(Q), a and b for the B c — > D* decay are listed in Table [J] 





F(0) 


a 


b 


V 


0.54 =f 0.018 


-1.28 


-0.230 


A! 


0.30 =F 0.017 


-0.13 


-0.180 


A 2 


0.36 =F 0.013 


-0.67 


-0.066 


« (A 3 - A ) 


-0.57 =F 0.040 


-1.11 


-0.140 


Ti 


0.31 =F 0.017 


-1.28 


-0.230 


Ti 


0.33 =F 0.016 


-0.10 


-0.097 


Ts 


0.29 =F 0.034 


-0.91 


0.007 



TABLE I. B c meson decay form factors in the three point QCD sum rules. 

The differential branching ratio is calculated without resonance contributions, including uncertainty on 
form factors, and with resonance contributions, and s dependence for 1/R — 200, 350, 500 GeV is presented 
in Figs. [5]and [21 respectively. The change in differential decay rate and difference between the SM results 
and new effects can be noticed in the figures. The maximum effect is around s = 0.25 ±0.05 (0.37±0.02) for 
H (r) in Fig. [2l In spite of the hadronic uncertainty, for 1/R — 200 GeV and 350 GeV, studying differential 
decay rate can be a suitable tool for studying the effect of extra dimension. 

Supplementary of these, 1/R dependence of differential branching ratio at s= 0.18 (0.4) for // (r) is 
plotted in Fig. 2J Considering any given bounds on the compactification factor the effect of universal extra 
dimension can be seen clearly for low values of R, with and without resonance contributions. On the other 
hand, 1/R > 600 GeV the contribution varies between ~ 5 — 8% more than the SM results. 

To obtain the branching ratio, we integrate Eq. (TT3]) in the allowed physical region. While taking 
the long-distance contributions into account we introduce some cuts around J/tp and ?/>(2s) resonances to 
minimize the hadronic uncertainties. The integration region for q 2 is divided into three parts for /x as 4m^ < 
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q 2 < (m J/4 , - 0.02) 2 , (m Jtf + 0.02) 2 < q 2 < (m M2s) - 0.02) 2 and (m^ (2s) + 0.02) 2 < q 2 < (m Bc - m^) 2 
and for r we have 4m 2 < q 2 < (m^^ 2s ) — 0.02) 2 (m^^s) +0.02) 2 < q 2 < (ttib c — m£>*) 2 , the same as in [48j . 

The results of branching ratio in the SM with resonance contributions and uncertainty on form factors, 
we obtain 



Br(B c -> £>>+ 
Br(B c -> D*t + t 



D* s n + f r) = 2.131°^ x 10 



j+0.27 
1.25 
h0.15 
-0.14 



■7 

1.4!) ur 8 . (15) 

Observing the contribution of the ACD, the 1/R dependent branching ratios, including resonance contri- 
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FIG. 5. (color online) The dependence of branching ratio on 1/R, including resonance contributions and uncertainty 
on form factors. 

butions and uncertainty on form factors, are given in Fig. [SJ Comparing the SM results and our theoretical 
predictions on the branching ratio for both decay channels, the lower bound for 1/R is found approximately 
250 GeV which is consistent with the previously mentioned results. 

As 1/R increases, the branching ratios approach to their SM values. For 1 /R > 550 GeV in both channels, 
they become less than 5% greater than that of the SM values. Between 1/R = 250 — 350 GeV the ratio 
is (2.66 - 2.40)lois x 10 ~ 7 for ^ (1-75 - l-61)±o:iJ x for T decay. Comparing these with the SM 

results, the differences worth to study and can be considered as a signal of new physics and an evidence of 
existence of extra dimension. 



IV. FORWARD-BACKWARD ASYMMETRY 



Another efficient tool for establishing new physics is the study of forward-backward asymmetry. The 
position of zero value of Afb is very sensitive to the new physics. The normalized differential form is 
defined for final state leptons as 

a / \ JO dsdz J —1 dsdz 

Afb{s) - r^dz+r ^dz 

JO dsdz J —1 dsdz 

(16) 

where z = cos and is the angle between the directions of £~ and B c in the rest frame of the lepton pair. 
In the case of B c -> D*£ + £~, we get 

G 2 a 2 m Bc |2 8m% r VXvs(Re[BE*} + Re[AF*]) 

A FB= 12 5 \VtbV ts 



2 12 vr 5 ' Ml dT/ds 
i% c V\vs{Re[BE*] + Re[AF*}) 



(17) 



Using above equation, we present the variation of lepton forward-backward asymmetry with s including 
uncertainty on form factors in Fig. [5] As 1/R gets smaller, there appears considerable difference between 
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the SM and the ACD results for s 5 0.16 in \i and 0.33 ;$ s < 0.43 in r decays. Considering the resonance 
contributions, the results are given in Fig. one can recognize a similar situation for s < 0.23 and 
0.32 < s < 0.44, respectively. 
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FIG. 6. (color online) The lepton forward-backward asymmetry including uncertainty on form factors. 
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FIG. 7. (color online) The lepton forward-backward asymmetry including resonance contributions 



To understand the dependence of Afb on 1/R for both lepton channels better, we perform calculation 
at s = 0.05 (0.4) for /Lt(r) and present the results in Fig. [5] In the fi channel, UED contribution on 
Afb gets important between 1/R = 200 — 600 GeV, while in r decay the contribution is insignificant for 
1/R > 400 GeF. 

The position of the zero of forward-backward asymmetry, Sq, is determined numerically and the results 
are presented in Fig. [5] Both plots for B c — > D*/j, + /i~is for the zero point in the s < 0.1 region; the lower 
(upper) one is for the resonance (non resonance) case, while the zero point for B c — > D*t + t~is because 
of resonance contributions. In the SM, resonance shifts the zero point of the asymmetry, so = 0.079, to 
a lower value, s = 0.068, in B c — > i.e., further corrections could shift s to smaller values [12| . 

As 1/R — > 200 GeV the sq approaches low values for both decay channels. In the 1/R = 250 — 350 GeV 
region, sq varies between (0.058 — 0.068) without resonance contributions and (0.051 — 0.058) with resonance 
contributions. The so shift is ~ 5% of the SM value for 1/R > 600 6*61^. The variation of so for B c — > 
D*t + t~is negligible. 



V. LEPTON POLARIZATION ASYMMETRIES 



We will discuss the possible effects of the ACD model in lepton polarization, a way of searching new 
physics. Using the convention followed by previous works (46|-(47j. in the rest frame of £~ we define the 
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FIG. 8. (color online) The dependence of lepton forward-backward asymmetry on 1/7? at s = 0.05 for /i and s = 0.4 
for r. 
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FIG. 9. (color online) The variation of zero position of lepton forward-backward asymmetry with 1/R. 



orthogonal unit vectors S i , for the polarization of the lepton along the longitudinal, transverse and normal 
directions as 



S~ = (0,e L )= 10,^ 
S- = (0,e N )= (0 



\Pe\, 

PD t x p e 



\Pd* xpe\J 

S T = (0,e T ) = {0,e N xe L ), (18) 

where pe and po- are the three momenta of £~ and D* meson in the center of mass (CM) frame of £ + £~ 
system, respectively. The longitudinal unit vector is boosted by Lorentz transformation, 

&=( M ^). (19) 

while vectors of perpendicular directions remain unchanged under the Lorentz boost. 

The differential decay rate of B c — > D*£ + £~ for any spin direction ri~ of the £~ can be written in the 
following form 



dr(rT) _ 1 fdT\ 



ds 2 \ ds J 

Here, (dT/ds)o corresponds to the unpolarized decay rate, whose explicit form is given in Eqn. ([1 



(20) 
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The polarizations P L , P T and P N in Eq. ([20)) are defined by the equation 



3r( n = )-3r( n 



-(n =e 2 ) + -(„ =-e s: ) 



for i = L, N, T. Here, P L and P T represent the longitudinal and transversal asymmetries, respectively, 
of the charged lepton t~ in the decay plane, and P^ is the normal component to both of them. 
The explicit form of longitudinal polarization for £~ is 



P7 = 



1 



-4m B v 



8mt sXRelAE*} + -<12rs + X)Re[BF*] 
r 



3A is 

—Xm% (1-r-s) \lle\BG*] + RelCF*]] + -X 2 m% RelCG*} 
r L J r 



(21) 
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FIG. 10. (color online) The dependence of longitudinal polarization on s without resonance contributions using 
central values of form factors. 
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FIG. 11. (color online) The dependence of longitudinal polarization on s with resonance contributions using central 
values of form factors. 



Similarly, the transversal polarization is given by 
1 



Ad? 



- 8m| Re[AB*} + ^ — - — '^-Re[BF*} - H^i? e [CF*] 
rs rs 



2 4 

m m 
— ^(1 - r)(l - r - s)Re[BG*] + — ^A(l - r)Re[CG*} 
rs rs 
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FIG. 12. (color online) The dependence of longitudinal polarization on 1/R including uncertainty on form factors 
and resonance contributions. 



B c 



(1 - r - s)Re[BH*] 



-Re[CH*} 



(22) 



and the normal polarization by 
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FIG. 13. (color online) The dependence of transversal polarization on s without resonance contributions using 
central values of form factors. 



m B minvv s 



4Im[BE*} - AIm[AF*] + -(1 - r - s)Im[FH*] 



--(1 + 3r - s)Im[FG*] - -m| \Im\GH*} 
r r 



(23) 



We eliminate the dependence of the lepton polarizations on s in order to clarify dependence on 1/R, by 
considering the averaged forms over the allowed kinematical region. The averaged lepton polarizations are 
defined by 



(Pi) 



Pi—ds 

(2m,/m Bc ) ! " s 



(l-m D »/m Bt .) ^ 



(24) 



(2me/m Bc ) 2 



ds 



ds 



The dependence of longitudinal polarizations on s with and without resonance contributions are given in 
Figs. [TU]and[TTl respectively. For high values of s as 1/R approaches 200 GeV the deviation from the SM 
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FIG. 14. (color online) The dependence of transversal polarization on s with resonance contributions using central 
values of form factors. 
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FIG. 15. (color online) The dependence of transversal polarization on 1/R, including uncertainty on form factors. 



results get greater for r in both resonance and non-resonance cases, while for /i channel this effect can be 
seen clearly for all s values when resonance contributions are not added; including resonance contributions, 
around the peaks this effect seems to be suppressed and only for low values of s we can mention a deviation. 
Eliminating the dependence of polarization on s, we get variation of longitudinal polarization with respect 
to 1/R, given by Fig. [T^l For 1/R> 500 GeV, the difference becomes less important for both channels. 
The SM longitudinal polarization, P L = -0.599, develops into -0.670 (-0.646) for 1/R > 250(350) GeV 
for /i. A similar aspect can also be noticed for r. That is, Pl — —0.321 SM value vary to —0.366 (—0.347) 
for 1/R > 250(350). 

The dependence of transversal polarization on s with and without resonance contributions are given in 
Figs. fl~3l and fT4l respectively. The UED effect is unimportant in both decay channels. In view of 1/R 
dependency, given by Fig. Qj)l no difference is observed for r decay. Up to 1 /R = 600 GeV the change is 
sizeable for /i channel. In particular, between 1/R. = 250 — 350 GeV the difference might be checked for a 
signal of new physics. 

We have plotted the variation of normal polarizations on s with and without resonance contributions in 
Figs. fT6land[T7l respectively and on 1/R in Fig. [18] The SM value itself for /i is tiny and as can be seen 
from the figures the effect of UED on normal polarization in this channel is irrelevant. Additionally, the 
relatively greater value of normal polarization in the SM for r differs slightly. 
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FIG. 16. (color online) The dependence of normal polarization on s with resonance contributions using central 
values of form factors. 
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FIG. 17. (color online) The dependence of normal polarization on s without resonance contributions using central 
values of form factors. 



VI. CONCLUSION 



In this work, we discussed the B c — > D*£ + £~decay for fj, and r as final state leptons in the SM and the 
ACD model. We used form factors calculated in QCD sum rules and throughout the work, we reflected the 
errors on form factors on calculations and demonstrate the results in possible plotting. 

Comparing the SM results and our theoretical predictions on the branching ratio for both decay channels, 
we obtain the lower bound as 1/R ~ 250 GeV. Although this is consistent with the previously mentioned 
results, a detailed analysis, particularly with the data supplied by experiments, is necessary to put a precise 
bound on the compactification scale. 

As an overall result, we can conclude that, as stated previous works in literature, as 1/R — > 200 GeV 
the physical values differ from the SM results. Up to a few hundreds GeV above the considered bounds, 
1/R > 250 GeV or 1/R > 350) GeV, it is possible to see the effects of UED. 

Taking the differential branching ratio into consideration, for small values of 1/R there comes out essential 
difference comparing with the SM results. 

Difference between the SM and the ACD results in the forward-backward asymmetry of final state leptons, 
particularly in the specified region, the obtained result is essential. In addition, the position of the zero of 
forward-backward asymmetry, which is sensitive in searching new physics, can be a useful tool to check the 
UED contributions. 

Polarization of the leptons have been studied comprehensively and we found that transversal and normal 
polarizations are not sensitive to the extra dimension, only dependence of transversal (normal) polarization 
on 1/R for fi (t) decay channel for low values of 1/R might be useful. However, studying longitudinal 
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FIG. 18. (color online) The dependence of normal polarization on 1/7? including uncertainty on form factors and 
resonance contributions. 



polarization for both leptons up to 1/R = 600 GeV will be a powerful tool establishing new physics effects. 

Under the discussion throughout this work, the sizable discrepancies between the ACD model and the 
SM predictions at lower values of the compactification scale can be considered the indications of new 
physics and should be searched in the experiments. 
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Appendix A: Wilson Coefficients in the ACD Model 

In the ACD model, the new physics contributions appear by modifying available Wilson coefficients in the 
SM. The modified Wilson coefficients are calculated in [13]- 18 1 and can be expressed in terms of F(xt, 1/-R) 
which generalize the corresponding SM functions Fo(x t ) according to 

oo 

) (Al) 



n=l 

where x t = ml/m^r, x n = m^jrtv^ and m„ = n/R. 

Instead of C7, an effective, normalization scheme independent, coefficient C7 in the leading logarithmic 
approximation is defined as 

C 7 e// ( M6 ,l/i?) =v 16/23 C 7 (^ w ,l/R) 



3 



| (r? 14/23 „ 7? 16/23 )c , 8(AiWj 1/R) + C2{ ^ Wi l/R) J2 krf* (A2) 



with 77 = aA f w } and 

' rK„f lit, 1 



a s (m z ) 



where in fifth dimension a s (mz) — 0.118 and /3q = 23/3. 
The coefficients a, and h; are 



The functions in (IA2I) arc 



Here, D'(xt, 1/-R) and E'(xt, 1/-R) are defined by using (|A1[) with the following functions 

3x? + 5x 2 -7x t ) x 2 (2-3x t ) 
12(1 -x t ) 3 + 2(1 -a:*) 



In ■ 

6(x t - l) 2 1 + x 



(2 - 7x„ + 3x 2 



(A3) 



/ 14 16 6 12 \ 
Gl = V23 ' 23 ' 23 ' ~ 23 ' °- 4 ° 86 ' -°- 4230 ' -°- 8994 ' °- 1456 ) 

hi = ^2.2996, -1.088, ~, -0.6494, -0.0380, -0.0186, -0.0057). (A4) 



C 2 (mw) = 1, C 7 (/i W r,l/7i) = -i£) / (^,l/ie), C 8 (ti w ,l/R) = ~E'{x t ,l/R). (A5) 



^ (a:t) = _12±i_^±i_^Z + ^-^J lnxt (A6) 



g °^=- ^^ + 2ra fa " (A7) 

, _ x t (-37 + 44x t + 17x 2 + 6x 2 (10 - 9x t + 3x 2 ) - 3x„(21 - 54x t + 17x 2 )) 

36(x t - \y 

(-2 + x„ + 3x t )(x t + 3x 2 + x 2 (3 + x t ) - x„(l + (-10 + x t )x t )) x„ + x t 



Z^r -~" In — — (A8) 

6 l + x„ y ' 

, x t (-17-8x t + x 2 -3x„(21-6xt+x 2 )-6x 2 (10-9x t + 3x 2 )) 

E ^ Xn) = i^—ir 

(1 + x„)(x t + 3x 2 + x 2 (3 + x t ) - x„(l + (-10 + x t )x t )) x„+x t 

H ; n m 

2(x 4 -l) 4 l + x„ 

-ix„(l + x„)(-l + 3x n )ln-^^. (A9) 

2 1 + X n 
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Following [TtJ or directly from [l2| one gets the expressions for the sum over n as 
^ ^(-37 + ^(44+17x0) 



ttM w R 



72( Xt - iy 



V v/2+7 f +%5/ W M,,iva 



(-2 + 3x t )x t (l + 3x t ) 
6(x t - l) 4 



J(R,-i/2) 



6(x t 


-1)4 




1 


6(x t 


-1)4 


(3- 


hX t ) 



+ 3x t ) - (-2 + 3a? t )(l + (-10 + x t )x t )]J(R, 1/2) 
[(-2 + 3as t )(3 + xt) - (1 + (-10 + x t )x t )]J(R, 3/2) 

J(R, 5/2) 



and 



6(x t - l) 4 
. x t (-17 + (-8 + i t )a*) 



24(x t - l) 3 



irM w R 



dy [y 1 ' 2 + 2z// 2 - 3y 5 / 2 ) coth(7rM w i2^) 



[xt - l) 4 

— j[ar t (l + 3xt) - (1 + (-10 + x t )x t )]J{R, 1/2) 



1 

(3 + it) 
(** " I) 4 



[(3 + x t ) - (1 + (-10 + xt)x t )]J{R, 3/2) 
J(#, 5/2) 



where 



J{R, a) = ( dy y a [coth(?r M w R^/y) - x\ +a cofa{nm t Ry/y)] • 
Jo 



The Wilson coefficient Cg in the ACD model and the NDR scheme is 

,NDR , Y{x u l/R) 



C 9 (n,l/R) = P<?»" + 



sin 2 8 



w 



4Z(x t ,l/R) + P E E(x t ,l/R) 



(A10) 



(All) 



(A12) 



(A13) 



where Pq DR = 2.6 ± 0.25 and Pe is numerically negligible. The functions Y(x t , 1/R) and Z(x t , 1/R) are 
defined as 



Y(x u 1/R) = Y (x t ) + C n (x t ,x n ) 

n=l 

oo 

Z(x t , l/R) = Z Q (x t ) + ^ C n (xt,x n ) 



with 



Y (x t ) = | 



3x t 



x t -A ^ 

x t - 1 (x t - I) 



-lnx t 



(A14) 



(A15) 



(A16) 
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Zo(x t ) 



18xj - 163a; 3 + 259x 2 t - 108x t 
144(x t - l) 3 



32xf - 38:r? - 15x2 + 18x t 1 



72{x t - l) 4 



9 



lnx t (A17) 



.T, 



0* " I) 2 



— 8xt + 7 + (3 + 3x t 4- 7x n — xtx n )ln 



Xt ~j~ Xyi 

1 + x n 



and 



£ ^o**' *») = i6 ( ( i t _?) - i6(fr^ 2 [3(i + 2:0 ~ i/2) + (a;t ■ ?) w i,2)] - 



The (U independent Cio is given by 



where Y(x t ,l/R) is defined in (IA14|) . 



Clo(W . JWffl 



(A18) 



(A19) 



(A20) 
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